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ABSTRACT 



A three-dimensional, primitive equation model by Anthes 
and Warner (1978) was used to predict five typhoon cases. 

The fxve-layer model on a 40X40 staggered grid with 120 km 
resolution included a moisture cycle, sensible and latent 
heat flux at the earth's surface, and a bulk parameteriza- 
tion of the planetary boundary layer. The model is initialized 
using operational wind fields and two forms of a wind bogus. 

Due to the lack of a representative moisture analysis 
in the vicinity of tropical cyclones, simulated moisture 
fields were used to initialize the model. Initial experiments 
conducted with these fields produced widespread convection 
and heating which developed circulations in areas well removed 
from the actual storm. The associated raodif icatio.ns to the 
steering flow, and the overly intense storm circulations 
resulted in premature recurvature. Use of a second moisture 
bogus, which provided less available moisture (especially at 
upper levels) , reduced the amount of convection over the 
entire grid and the effect on the steering current. 

A second wind bogus based on a scheme described by Hover- 
male (1976) was also tested. This technique resulted in a 
smaller and less intense initial st.orm, which also resulted 
in considerably less convection. The tracks forecast by 
the different wind and moisture fields were compared with 
forecasts of the 60 km resolution Madala and Hodur model 
using the same initial data. 
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I. INTRODUCTION 



The major objective of this set of experiments was 
to test a general mesoscale model using operationally- 
analyzed data to determine effectiveness in forecasting 
tropical cyclone development and movement compared to some 
other tropical cyclone models currently being tested. The 
Penn State model was selected because it contained boundary 
layer physics, a convective moisture cycle, and other 
options which could be easily included in (or excluded 
from) any forecast. It was then possible to determine the 
effects and relative importance of the various parameteriza- 
tions by selectively incorporating different combinations. 
This model was also chosen for its general code which allowed 
changes to be made easily to the resolution and size of 
the vertical and horizontal domains. These changes were 
necessitated by time and computer availability constraints. 

During the initial forecasts, it was discovered that 
the initial wind and moisture conditions played a very 
important role in determining the development and movement 
of tropical cyclones. An investigation was then conducted 
to determine the extent of the effects of the initial wind 
and moisture conditions on the model forecasts. This paper 
presents the results of these experiments and compares fore- 
casts made using the Penn State model with forecasts made 
using two other tropical cyclone models. 
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II. DESCRIPTION OF THE MODEL 



The experiments used a general, three-dimensional, 
primitive equation model designed to forecast flows with 
characteristic wavelengths of 10-2500 km. Provisions 
of the model included a variable terrain, a moisture cycle, 
latent heating, sensible heat flux at the earth's surface 
and boundary layer physics. The moisture cycle modeled 
both stable (nonconvective) precipitation with grid-scale 
relative humidities in excess of 100%, and unstable (convec- 
tive) precipitation with grid-scale relative humidities 
less than 100%. Details, such as the finite difference 
equations and staggered grid lattice, were retained intact 
as described by Anthes and Warner (19 78) . 

The particular version of the model used in this set 
of experiments was similar to one described by Anthes 
(1978). Model inputs consisted of pressure weighted wind 
components (p*u,P*v) , relative humidity and temperature at 
five levels, surface pressure, sea surface temperature and 
terrain. Predicted parameters were wind components (u,v) , 
surface pressure (p^) , relative humidity (moisture) , and 
potential temperature (9). Anthes ' s (1978) version was 
reduced from 7 levels to 6 levels at which a was defined 
(a = 0.0, 0.3, 0.5, 0.7, 0.9 and 1.0) with all other variables 
(V,T,<t>,a), and q) defined halfway between these levels at 
a = 0.15,0.4,0.6,0.8 and 0.95. Sigma level 0 . 0 corresponded 
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to a pressure of 0 millibars. To reduce computer time, 
a 40x40 grid in Mercator coordinates with 120 km resolution 
was selected rather than the 50x50 grid with 60 km resolution 
used by Anthes (1978) . Other simplifications were made to 
the model for economic reasons as described by Anthes (1978) . 
A spatially and temporally constant vertical distribution 
of latent heating associated with cumulus convection was 
specified in lieu of a one-dimensional cloud model (Anthes , 

1977) , and a simplified version of Deardorff's (1972) bulk 
parameterization of the planetary boundary layer (PBL) was 
substituted for the high resolution PBL (Anthes and Warner, 

1978) . This scheme assumed that the surface wind, Iv I, 
was equal to the wind at the lowest level (o = 0.95). 

All experiments described here used an ocean terrain with 

-3 

no land, and as such, a constant value of 1.5x10 was 
used for the drag coefficient, 

All layers were checked for supersaturation and 6 and 
q were adjusted to maintain relative humidities less than 
100% by condensing a fraction of the moisture excess and 
adding the associated latent heat to the temperature ten- 
dency. The original model assumed a constant fraction of 
.31 for all levels. This factor was modified using an 
approximation suggested by Haltiner (personal communication) . 
A new constant was calculated for each level using the 
equation 



Aq 



q3 - q 



(l + L^q^/Cp T^I 



( 1 ) 
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where typical values for q and T were obtained from 
Jordan's (1958) mean September tropical sounding and 
Sheets' (1969) mean hurricane sounding. The modified 
values which were equal to the reciprocal of the denomina- 
tor in equation (1) are listed in Table 1. 

Initial conditions were derived from the Fleet Numerical 
Weather Central (FNWC) Global Band Analysis. Initiali- 
zation was accomplished using a scheme developed by Hodur 
(personal communication) for use with the Madala and Hodur 
(1977) split semi-implicit tropical cyclone model. This 
scheme used a reverse balancing technique to derive interior 
geopotential values on sigma levels from the nondivergent 
wind components using a nonlinear balance equation. Balancing 
was conducted over a 51x51 grid with 120 km resolution, 
from which the 40x40 grid was interpolated for input to 
the model. 



Table 1. Fractions used for condensing excess 
moisture to obtain 100% relative 
humidity. Values for T and q were 
taken from representative tropical 
soundings . Values of q^ were calcu- 
lated using the Clausius-Clapeyron 
equation. 



Level (a) 


T (°K) 


0 . 15 


206.76 


0 . 40 


258.26 


0.60 


275.66 


0 . 80 


286.36 


0.95 


295.76 



qg (gkg ^) 


Fraction 


0.0433 


0.986 


3.0399 


0.619 


7.6894 


0.422 


12.0868 


0.335 


18.8008 


0.256 
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Iriitial wind conditions were composed of both analyzed 
and idealized values . Following Hodur (personal communica- 
tion) , a symmetrical wind bogus (referred to as V— bogus 1) 
was obtained within a 360 km radius of the storm center by 
3-dding idealized values to the vorticity field (which was 
derived from the analyzed winds) prior to calculation of 
the stream function and nondivergent winds. Bogus vorticity 
values were calculated for the surface to 700 mb from 



? (k) 



2V 



r 



max 

o 



( 2 ) 



where was the magnitude of the wind at the radius 

of maximum wind (^q) • The vorticity maximum at 100 mb was 
set equal to 1.0x10 ^ s with the vorticity at inter- 
mediate levels varied linearly between 700 mb and 100 mb. 
Horizontally, the vorticity was decreased in accordance 
with 



?^(k) = C_ 3 .^(k) , 0 < r < 360 

(3) 

where r was the radial distance (km) from the storm center. 
Values for were taken from Typhoon Reports and typhoon 

warning messages and the radius of maximum wind (RMW) , r^, 
was set equal to 60 km. 

Due to the lack of representative moisture analyses in 
the vicinity of tropical cyclones, idealized moisture fields 
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(referred to as Rh-bogus 1) were used to initialize the 
model. These fields were produced by specifying a uniform 
field of large scale relative humidities which decreased 
with height. Further bogusing of the storm was accomplished 
by centering a five grid point square of higher relative 
humidities, which also decreased with height, over the 
storm. In this manner, the storm started with more mois- 
ture than the surrounding environment. Large scale and 
storm relative humidities are given in Table 2. 



Table 2. Relative 

Level (a) 

0.15 

0.40 

0.60 

0.80 

0.95 

Relative 
Level (a) 



0.15 

0.40 

0.60 

0.80 

0.95 



humidities used in Rh- 

Storm 

.80 

.80 

.90 

.90 

.95 

humidities used in Rh 



Storm 



low 


high 


. 30 


.40 


.53 


.60 


.76 


.87 


.77 


.87 


. 80 


.90 



bogus 1 and Rh-bogus 2 

Large Scale 

.50 

.60 

.70 

.80 

.85 

bogus 2 . 

Large Scale 



.20 
.46 
.51 
.69 
. 80 
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III. EXPERIMENTS 



The basic experiment consisted of 48-hour forecasts 
for the following five cases using all of the model's 
options (with the exception of terrain) ; 

I. Typhoon OLIVE 78042112 (year , month, day , 12 GMT) 

II. Typhoon OLIVE 78042400 

III. Typhoon KIM 77110912 

IV. Typhoon PAMELA 76051900 

V. Typhoon PAMELA 76052200. 

Tracks of these storms were then compared to the tracks 
obtained from forecasts made with the Madala and Hodur 
(1977) split semi-implicit (SSI) model, the FNWC Opera- 
tional Tropical Cyclone (Channel) model (Mihok and Kinsman, 
1977), and the best track analyses. Appendix A contains a 
comparison of these two models with the Penn State model. 
Further experiments were then conducted to test the model's 
sensitivity to the initial wind and moisture conditions and 
to the convective parameterization. 

A. RESULTS OF THE FIRST COMMON EXPERIMENT 

All five cases experienced strong recurvature. See 
Figs. 17 through 21 for positions relative to the best track. 
Examination of the wind and moisture fields revealed several 
characteristics which were common to all cases. Every storm 
developed strong frictional convergence with large cross- 
isobaric flows at low levels, even in areas of high wind 
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velocities. Although the maxiim.im storm intensities varied, 
all five cases developed a southerly flow at middle levels 
just to the east of the storm. An example of this circu- 
lation, which was especially prevalent at 600 mb for Typhoon 
KIM, can be seen in Fig. 1. A tendency to produce a large 
anticyclonic circulation at 150 mb with strong westerlies 
located along the northern branch of the gyre also existed 
in all cases. An example of this circulation is given in 
Fig. 2. This circulation, combined with the strong middle 
level southerly flow, probably contributed significantly 
to recurving the storms. 

Another dominant characteristic common to all cases, 
which contributed to the recurvature effect, was the produc- 
tion of widespread convection, with particularly strong 
convective activity occurring after 12 to 24 hours in the 
northeast corner and along the northern boundary. All 
cases experienced excessive vertical velocities accompanied 
by strong latent heating and precipitation (both convective 
and nonconvective) in these areas, which resulted in the 
formation of localized areas of lower pressures. Continued 
development of a trough in the northeast corner distorted 
the entire surface pressure pattern to such an extent that 
the storm was steered toward the corner, producing even 
greater development in this region and more rapid recurva- 
ture of the storm. An example of this sequence can be 
seen in Figs. 3a through 3d. Forecasts for both cases of 
Typhoon OLIVE also developed severe convection in the 
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southwest corner and an intense circulation was produced 
between 24 and 36 hours. The first OLIVE case also 
developed an additional circulation which had existed 
initially to the west of the storm. An example of the 
areas which underwent development due to convective activity 
can be seen in Figs. 4 through 6. The development of these 
circulations near the boundaries prompted an investigation 
into the moisture fields, latent heat parameterization 
scheme and the boundary conditions. These topics are 
examined in the following sections. The general procedure 
was to select a typical case and examine successively the 
effect of each modification on that storm. After adopting 
the modifications based on the one case, the entire set of 
five cases was re-run. 

B. MOISTURE EXPERIMENTS 

A new moisture initialization was introduced in an attempt 
to reduce the amount of convection and high relative hiimidi- 
ties at upper levels in regions not associated with the 
storm. The second technique (referred to as Rh-bogus 2) 
also consisted of an area of higher relative hximidities 
superimposed on a uniform field of lower values. However, 
the second moisture bogus covered an area with a radius of 
3 grid points (360 km) , which was the same area covered 
by the wind bogus. The large scale values, which were 
based on Jordan's (1958) mean September tropical sounding, 
varied only with height, while the storm relative humidities. 



20 



taken from Sheet's (1969) mean hurricane sounding, depended 
on the storm surface pressure as well as height. Values 
used in this scheme are listed in Table 2. The high range 
of relative humidities was used for storm surface pressures 
less than 99 5 m±) and the low set for pressures greater than 
1014 mb, with a linear interpolation between those sets 
for intermediate storm surface pressures. Note that the 
major difference in Rh-bogus 2 was the reduction of all 
relative humidities, especially at upper levels. Additional 
48-hour forecasts were then conducted for both cases of 
OLIVE using the new moisture bogus. OLIVE was chosen 
because it had the most widespread convection of all the 
cases . 

C. RESULTS OF RH-BOGUS 2 

The results of the low moisture bogus are typified 
in Figs. 7 through 9. Convective activity was reduced over 
the entire grid with the greatest effect noticeable in the 
reduced development of those areas not associated with the 
storm. The total convective rain produced over the grid 
was reduced from 382 cm to 198 cm, while the total non- 
convective precipitation was reduced from 26 cm to 7 cm 
in the first 12 hours. More significant, however, was the 
reduction of the maximum convective precipitation associated 
with the second cyclone (located to the west of OLIVE) 
from 12.5 cm to 1.4 cm. The maximum convective precipitation 
associated with OLIVE was only reduced from 5.9 cm to 4.5 cm. 
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The associated convective heating rate at 12 hours was 
lowered from 280 °C/day to 60 °C/day over the second 
cyclone, while the heating over OLIVE was reduced from 
160 °C/day to 110 °C/day. Similar results were obtained 
in the other areas, particularly in the northeast corner, 
although the development of the trough in this region was 
not prevented. 

The net effects of the low moisture bogus were filling 
the storm by about 2 mb and decreasing the winds by less 
than 10%, while significantly reducing the deepening over 
the large scale areas. The general characteristics of the 
flow were not altered, nor were the forecast tracks signi- 
ficantly altered. 

D. NO HEAT EXPERIMENT 

A 24-hour forecast was then made using the high humidity 
bogus, but without including latent heating effects in the 
temperature prediction equation. Typhoon KIM was chosen for 
this experiment because it had the best predicted track of 
the five cases using the SSI model and had been used with 
and without the wind bogus (Hodur, personal communication). 
Deletion of latent heating effects (due both to convective 
and nonconvective processes) resulted in rapid filling of 
the storm and only marginal improvement of the storm track. 

E. CHANGES TO THE MODEL 

Before continuing with further experiments, several 
detailed changes were made to some of the model-defined 
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constants to make the 5~layer version more closely resemble 
the original 6-layer model described by Anthes (1978) . 

The vertical distribution of convective heating, N(a) , 
defined by Anthes (1977) , was altered in the lowest two 
layers such that the same amount of heating was provided, 
but over different layer thicknesses. Because this function 
was normalized, the changes produced in these two layers 
were distributed throughout the vertical. The corrected 
values of N(a), which provided the same distribution of 
convective heating described by Anthes (1978) in the 6- 
layer model, are compared in Table 3. A 24-hour forecast 
using the high moisture bogus and the revised heating 
distribution produced a 10% reduction in convection and 
its associated latent heating and precipitation. The srorm 
center filled 1 mb. 

Two other changes were made to the convective parameteri- 
zation. The vertical eddy fluxes, 6o)^ (q^-q) /5a , which 
determined 3q/3t in the cloud, were altered to make the 
vertical distribution of cloud moisture more representa- 
tive of the distribution outlined by Anthes (1977). A 
24-hour forecast produced less convection, precipitation, 
and heating, and increased the surface pressure of the storm 
by 1 mb. Convection over the entire grid was further 
reduced in a 12-hour forecast by increasing the critical 

value of the integrated water vapor convergence, M^, required 

— 7 — 6 

for convection to occur from 3.0x10 to 3.0x10 9 ^ 2 ^^ ^ 

where 
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Table 3 



. Revised values of the vertical 
distribution of latent heating, 
N(a). (New values are given 
for experiment designed to shift 
heating maximum to higher levels.) 



ivel (mb) 


6 -layer 


5-layer 


New 


150 


- 


. 375 


.7595 


300 


1.065 


- 




400 


- 


1.5 


1.5190 


460 


1.308 


- 




580 


1.356 


- 




600 


- 


1.625 


1.2658 


700 


1.162 


- 




800 


- 


1.25 


1.0126 


820 


.872 


- 




940 


.194 


- 




950 


— 


.125 


.1266 



24 



1 

g 



M, 



I V(p*Vq)6aj^ 
k=k. 



+ i(p*qa) 



The number of grid points experiencing convection for any 
one time step was reduced significantly, with the maximum 
number of grid points (at 12-hours) lowered from 797 to 
only 222. The total convective precipitation over the grid 
was reduced from 414 cm to 263 cm. The changes in convec- 
tive precipitation and heating produced by the revisions 
discussed here, can be seen in Figs. 10a through lOd. 
Although the motivation for the above changes was to reduce 
what appeared to be excessive latent heating in the model, 
note that the circulation in the northeast corner still 
experiences strong convergence and excessive precipitation. 

Two experiments designed to test the model's sensitivity 
to the convective parameterization were conducted using Rh- 
bogus 2 for initial conditions. A 24-hour forecast was 
made using a vertical distribution of latent heating which 
shifted the maximum heating to higher levels . This was 
done so that the model would develop more cyclonic circula- 
tion at upper levels and possibly produce a greater westward 
movement. The new values of N(c) are listed in Table 3. 
Additional upper level heating produced less precipitation 
in the northeast corner but failed to alter the forecast 
significantly in 24-hours. Another 24-hour forecast was 
made with the integrated water vapor convergence, M^, 
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calculated for only the two lowest layers of the model. 

The objective of this experiment was to remove the effects 
of the large upper level moisture convergence and, thereby, 
reduce the precipitation in the northeast corner while not 
seriously affecting the storm. Confinement of the moisture 
convergence calculation to low levels reduced convection and 
development of a trough in the northeast corner while 
intensifying convection associated with the storm after 
12 hours. However, after 24 hours, the trough in the 
northeast corner developed rapidly and produced results 
similar to those already encountered. These last two 
model revisions were not incorporated into the model for 
future experiments. 

Frictional effects were also modified for changes made 
in 60 and p* which had resulted from changing the pressure 
at the top of the model from 200 mb to 0 mb, and decreasing 
the thickness of the surface layer from 150 mb to 100 mb. 
These modifications resulted in a 4% decrease in the effec- 
tive drag coefficient. A 12-hour forecast made using the 
reduced friction deepened the storm surface pressure by 
1 mb and increased the amount of convection occurring over 
the storm, while decreasing the total number of grid points 
experiencing convection. 

F. HORIZONTAL DIFFUSION EXPERIMENT 

An attempt was made to reduce the strong horizontal 
shears which existed along the boundaries. An example of 
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this shear can be seen in Fig. 2 where the interior upper 

level circulation differs drastically from the conditions 

along the western and southern boundaries. Since this 

characteristic also occurred when the model was run with 

extratropical data (see Anthes , 1978, Figs. 23c through 

23e) , it was decided to increase the horizontal diffusion. 

The horizontal eddy viscosity, was multiplied by an 

amplification factor which increased from a value of 1.0 

near the center to a maximum near the boundaries , where 

4 2-1 

Kjj was restricted to be less than 30 x 10 m s . All future 

runs with the model were made using twice the calculated 

4 2 -1 

Kj. and a maximum value of 60x10 m s . The increased 
diffusion decreased the shears along the boundaries by less 
than 10% and therefore did not significantly reduce the 
convergence which resulted in convection along the northern 
and eastern boundaries. 

Another factor which may have contributed to over- 
development of the northeast corner was the lack of moisture 
diffusion. The winds, especially at upper levels, tended 
to advect moisture into that corner. Without diffusion, 
it may have been possible that moisture simply accumulated 
in the corner producing more convection. This effect can 
be seen by comparing Fig. 3c with Fig. 11 which is the 24- 
hour forecast without latent heating. In the latter, although 
convection occurred, the heating effects were not added to 
the tendency and the circulation near the corner did not 
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develop. As a result, the surface pressure pattern was 
not distorted. 

G. WIND BOGUS EXPERIMENTS 

As demonstrated in Fig. 3a, the initial wind bogus 
(V-bogus 1) produced a storm which was very large. Although 
bogused values were entered in the vorticity field within 
a radius of 3 grid points, the cyclonic circulation and 
the associated pressure fields derived from the balance 
equation extended out to a radius of 8-10 grid points. 

It was also desirable to include an inner region of cyclonic 
outflow at upper levels. A second wind bogus (V-bogus 2) 
was calculated using the same conditions from the surface 
to 700 mb as V-bogus 1, while at 400 mb, the vorticity was 
given by 



= C 



max 



180 ^^max 



f) 



r < 180 km 



(4) 



r = -f r = 180 km , (5) 

^r 

c = -f + ~ f 180 < r < 360 km, ( 6 ) 

^r 180 






r 



= 0 



r = 360 km, (7) 



where the vorticity maximum, 5 was computed in the same 

IT13.X 

manner as described earlier for V-bogus 1. At 100 mb, the 
vorticity was calculated using one-half the calculated 
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values for 400 mb. The vorticity at other levels (500, 
300,250,200, and 150 mb) was varied linearly between 700 
and 400 mb and between 400 and 100 mb. This bogus produced 
essentially the same results as V-bogus 1, probably because 
the circulation in the lowest levels of the storm was still 
very large. 

Therefore, a third wind bogus (V-bogus 3) was calculated 
in an attempt to reproduce a wind bogus described by 
Hovermale et al (1976). The tangential wind, Vg, was 
calculated from 




22 r 
60 



0 < r < 60 km, (8) 



and 





60 < r 



1 



(9) 



for the lowest level, while the upper level wind was 



derived using 




12 r 
60 



0 < r < 6 0 km 



( 10 ) 




12 [1 - (%^)] 60 < r < 660 km. 



( 11 ) 
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and 



Vg (m s 



Tori _i_ km ) 1 ^ ^ ^ 

■" 600 km ' ««° i 



( 12 ) 



where r extended out to a radius of 6.4 grid points (768 
km) . The vorticity was then calculated from 






!li + !i - i !l£ 

3r r r 90 



(13) 



where the last term was taken as zero for a symmetrical 
storm. Prior to substituting the bogused vorticity values, 
a simple Laplacian smoother was passed 30 times over the 
analyzed wind components in the same area covered by the 
bogus. Experiments conducted on both V-bogus 1 and V- 
bogus 3 with and without the smoother did not produce 
significantly different results for either bogus. 

The following experiments were then conducted using 
Typhoon KIM and Rh-bogus 2 ; 

I. V-bogus 1: 24-hours 

II. V-bogus 3: 48-hours 

III. V-bogus 1 with one-half the vorticity 
maximum at the surface as in case I: 

48-hours . 

Differences produced in the initial wind fields by the two 
bogusing techniques can be seen in Fig. 12a through 12d. 
Note that V-bogus 1 produced a very large and intense storm 
compared to the one produced by V-bogus 3. In particular. 
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V-bogus 1 produced large changes in the winds and pressures 
at grid points a considerable distance from the storm. 

Because the 40x40 grid was selected from the interior of 
a 51x51 grid, wind values were changed as much as 50% in 
magnitude and 180'’ in direction along the boundaries by 
V-bogus 1. Subsequent analysis fields, that were required 
for calculation of the time-varying boundary values, were 
not bogused. The differences in these two fields introduced 
accelerations greater than expected in a 12-hour period 
along the southern and eastern boundaries. These accelera- 
tions were not introduced by V-bogus 3. 

H. RESULTS OF WIND BOGUSING EXPERIMENTS 

After 12 hours, V-bogus 1 produced convection at 100-200 
grid points while V-bogus 3 produced convection at only 20-30 
grid points. Total convective precipitation over the grid 
was lowered from 150 cm to 96 cm and the total stable rain- 
fall was reduced by about 40%. V-bogus 3 dramatically reduced 
the amount of precipitation and heating in the northeast 
corner while concentrating more rain and heating over the 
storm. The effects of the smaller vortex can be seen in 
Figs. 13 through 15. Note the lack of development of a 
large trough in the northeast corner. The southerly flow 
at middle levels was also reduced. Reduced development of 
these two features resulted in the improved forecast track, 
shown in Fig. 16. 
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The 24-hour forecast developed more convection in the 
northeast corner, however the storm seemed to be suffi- 
ciently far from the corner to not be seriously affected 
by the weak troughing. A large anticyclonic circulation 
developed at 150 mb with strong southerly winds in the 
southwest quadrant. Apparently, these winds contributed 
to the northward movement of the storm. 

Reducing the vorticity of V-bogus 1 to half the 
original value {Experiment III) resulted in a smaller 
vortex which produced values of convective activity, precipi- 
tation and heating which were between those produced by 
V-bogus 3 and the high intensity V-bogus 1. The track fore- 
case using this bogus was also between the tracks produced 
by the other two boguses. This experiment demonstrated 
the effect that variations in the Coriolis parameter can 
have upon excessively large or intense vortices; that is, 
larger storms produce larger northward accelerations 
(Hovermale, et al . , 1977) . 



32 



IV. COMPARISON WITH OTHER MODELS 



The effects of V-bogus 1 were compared to the FNWC 
analyzed fields (no storm bogus) for KIM, OLIVE (78042400) 
and PAMELA (76052200) using the SSI model by Hodur (personal 
communication) . A 48-hour forecast was also made for 
PAMELA (76052200) using the high humidity moisture bogus 
(Rh-bogus 1) and V-bogus 1. A summary of the five cases 
available for comparison with experiments using the Penn 
State model is given in Table 4. 



Table 4. 



Summary of experiments using the 
various bogusing techniques. 



Storm SSI Penn State 



OLIVE 


(78042112) 


Analyzed winds, dry 


V-bogus 1/Rh-bogus 
V-bogus 3/Rh-bogus 


OLIVE 


(78042400) 


Analyzed winds, dry 
V-bogus 1, dry 


same 


as above 


KIM 




Analyzed winds, dry 
V-bogus 1, dry 


same 


as above 


PAMELA 


(76051900) 


Analyzed winds, dry 


same 


as above 


PAMELA 


(76052200) 


Analyzed winds, dry 
V-bogus 1, dry 
V-bogus 1, wet 


same 


as above 



Forecast tracks obtained from the common experiments are 
shown in Figs. 17 through 21. V-bogus 1 produced more 
recurvature in all cases tested with the SSI model, although 



1 

2 
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not as much as it did in the Penn State model. The SSI 
model did not forecast the troughing that the Penn State 
model did in the northeast corner, although an area of 
lower surface pressure did exist there. The 12-hour fore- 
cast using the SSI scheme for KIM and PAMELA developed 
strong southerly flows to the east of the storm at middle 
levels which were similar to those produced by the Penn 
State model using V-bogus 1. However, the SSI model also 
developed strong northerly flows to the west of both storms 
at middle and upper levels. These winds acted to reduce 
the northward movement of the storm. The lack of the 
northerly winds, coupled with the trough in the northeast 
corner and the strong upper level anticyclone, allowed the 
southerly winds to move the storm a greater distance to 
the north in the Penn State model. 

V-bogus 3 and Rh-bogus 2 decreased the upper level 
anticyclonic flow and strong southerly flow in all cases 
using the Penn State model. Although all cases continued 
to experience convection in the northeast corner and along 
the northern boundary, development of a strong trough was 
reduced and limited to forecast times in excess of 24 
hours. Use of V-bogus 3 and Rh-bogus 2 resulted in improved 
tracks for all cases using the Penn State model. Insertion 
of V-bogus 1 in the SSI model reduced the errors in storm 
track positions for OLIVE and PAMELA and increased the 
errors for KIM (Hodur, personal communication). Rh-bogus 1 
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improved the forecast positions in the SSI model for PAMELA 
slightly over those using V-bogus 1 without moisture effects. 

The Penn State model produced smaller errors than the 
SSI model in storm positions at 12 hours for most of the 
five cases. After 12 hours, the Penn State model using 
V-bogus 3 and Rh-bogus 2 produced smaller errors than the 
SSI using FNWC analyzed winds (no bogus) , but larger errors 
than the SSI with V-bogus 1. The exception was Typhoon 
KIM for which all versions of the SSI scheme produced 
smaller errors than the Penn State model after 12 hours. 
Position vector errors are summarized in Table 5. 

Finer resolution was probably responsible for the over- 
all better accuracy obtained with the SSI model. The 
improved resolution allowed the storm to interact with 
small-scale circulations. With less resolution in the 
Penn State model, the large vortex storm was subjected 
primarily to the large-scale steering currents, and 
particularly more steering by the westerlies at upper levels, 
and consequently, more recurvature. 

Both cases of PAMELA produced large errors in later 
forecasts due to insufficient recurvature when V-bogus 3 
was used for initialization. The lack of recurvature proba- 
bly resulted from the rapid filling and deterioration of 
storm structure which existed in all experiments using 
PAMELA. Observational data indicated that the maximum 
wind speed in the storm was much more intense (60 m s ^) , 
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Table 5. Summary of position vector errors (errors are in nautical miles) 
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yet V-bogus 3 was designed to produce maximum winds of 
22 ms Accordingly, V— bogus 3 was altered to include 

more intense wind speeds by changing equations (8) and (9) 
to 

Vg (m s~^) = 60 ^ 0 < r < 60 km (14) 

and 

Vg(m s"^) = 60 - 60 < r (15) 

for the same limits of r. This bogus (referred to as V- 
bogus 4) assumed the same maximum wind and radius of maximum 
wind as V-bogus 1, but produced a significantly smaller 
storm, as can be seen by comparing Fig. 22 and Fig. 23. 

As a result of a greater vorticity maximiam and the smoother, 
V-bogus 4 produced greater maximum winds (67 ms ^) than 
V-bogus 1 (57 m s ^) . However, the greater winds did not 
increase the size of the storm. Because the storm was 
still reasonably small, boundary values closely resembled 
those of V-bogus 3 and the original analyzed values, whereas 
V-bogus 1 had created significant changes in boundary values. 
Resultant vorticities produced by the different bogusing 
techniques are plotted in Fig. 24. Note the faster decline 
in vorticity with increasing radius produced by V-bogus 3 
(and the similar V-bogus 4) . 
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48-hour forecasts made using V-bogus 4 and Rh-bogus 2 
produced dramatic improvements in the track for Typhoon 
PAMELA (see Fig. 25). By 24 hours, V-bogus 4 produced 
300-400% more convective precipitation and heating over 
the storm, but increased convection in other areas by less 
than 10%. 
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V. CONCLUSIONS 



The primitive equation model tested appeared to be 
very sensitive to latent heating effects produced by 
convection. Since the convective parameterization was 
initiated by the vertically integrated moisture conver- 
gence / proper representation of initial wind and moisture 
conditions at all levels was important. A moisture bogus 
consisting of a uniform field of low relative humidities 
with an area of higher relative humidities centered over 
the storm enhanced development of the storm while suppressing 
convection and development in areas not associated with the 
storm. The model also demonstrated some sensitivity to 
small changes in the initial values used in the moisture 
bogus . 

Experiments conducted using the Penn State model and 
another primitive equation model in conjunction with a 
simple vorticity bogus (V-bogus 1) showed that forecast 
storm tracks could be improved over forecasts made using 
the FNWC analyzed wind fields. Premature recurvature pro- 
duced in the Penn State model forecasts was retarded by 
using a smaller, less intense vortex (V-bogus 3). Alterna- 
tively, more recurvature was induced by increasing the 
maximum wind used in V-bogus 3 . This bogus produced more 
intense winds than V-bogus 1 without introducing a large 
cyclonic area. The smaller more intense vortex improved 
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forecast tracks over V-bogus 1 and the less intense version 
of V-bogus 3 in the Penn State model forecasts. It was 
therefore possible to produce significant, and somewhat 
controlled, changes in the forecast tracks of individual 
storms by varying the size and intensity of the initial 
bogus . 

Future experiments should include increased resolution 
of the horizontal domain. The SSI model appeared to have 
an advantage in representing steering currents because of 
its finer resolution. The initial wind maximum introduced 
by the bogus could not be sustained on the 120 km grid, 
and marked reductions in intensities occurred for all 
boguses used. Given finer resolution, V-bogus 3 using a 
variable maximum wind based on observations of the storm 
should provide a well-structured storm which does not get 
too large. 

The development of strong shears and convergence which 
produced convection along the boundaries , indicated that 
modifications to the boundary conditions should also be con- 
sidered in future experiments. Changes to the model which 
could resolve these effects are the inclusion of moisture 
diffusion and the addition of large-scale changes close to 
the boundaries. One possibility is the use of a scheme 
described by Perkey and Kreitzberg (1976) where a fraction 
of the lairge-scale time change (based either on analyses or 
from forecasts produced by a large-scale model) is added to 
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the model predicted tendency depending on the proximity 
to the boundary. 
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Figure 1. Typhoon KIM 24-h forecast VT 77111012 
600 mb winds (m s~l) and relative 
humidity (coordinates are grid points 
vice latitude and longitude) 
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Figure 2. Typhoon KIM 48-h forecast VT 77111112 
150 mb winds (m s“^) and relative 
humidity (isoline spacing is 20%) 
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Figure 3a. Typhoon KIM initial field SLP-1000 mb 
and winds (m s“^) at a = .95 using 
V-bogus 1 (isoline spacing is 4 mb) 
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Figure 3b. Typhoon KIM 12-h forecast SLP-1000 mb 
and winds (m s~^) at a = .95 using 
V-bogus 1 and Rh-bogus 1 
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Figure 3c. 



Typhoon KIM 24-h forecast SLP-1000 mb 
and winds (m s~l) at a = .95 using 
V-bogus 1 and Rh-bogus 1 
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Figure 3d. Typhoon KIM 36-h forecast SLP-1000 mb 
and winds (m s“^) at a = .95 using 
V-bogus 1 and Rh-bogus 1 
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Figure 4. Typhoon OLIVE 36-h forecast VT 78042300 
SLP-1000 mb and winds (m s~l) at a = .95 
using V-bogus 1 and Rh-bogus 1 
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Figure 5. Typhoon OLIVE 36-h foreoast VT 78042300 
accuraulated convective precipitation 
(cm) using V-bogus 1 and Rh-bogus 1. 
Center value of Olive is 20 cm (isoline 
spacing is 1 cm) 
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Figure 6. Typhoon OLIVE 36-h forecast VT 78042300 
convective heating (°C/day) using 
V-bogus 1 and Rh-bogus 1 (isoline 
spacing is 25 °C/day) 
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Figure 7. Typhoon OLIVE 36-h forecast VT 78042300 
SLP-1000 mb and winds (m s“^) at 
a = .95 using V-bogus 1 and Rh-bogus 2 
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Figure 8. Typhoon OLIVE 36-h forecast VT 78042300 

accumulated convective precipitation (cm) 
using V-bogus 1 and Rh-bogus 2 
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Figure 9. Typhoon OLIVE 36— h forecast VT 78042300 
convective heating (°C/day) using 
V-bogus 1 and Rh-bogus 2 
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Figure 10a. Typhoon KIM 12-h forecast VT 77111000 

accumulated convective precipitation (cm) 
using V-bogus 1 and Rh-bogus 1 before 
making changes to model constants (see text) 
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Figure 10b. Typhoon KIM 12-h forecast VT 77111000 

accumulated convective precipitation (cm) 
using V-bogus 1 and Rh-bogus 1 a'fter 
changes to model constants (see text) 
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Figure 10c. Typhoon KIM 12-h forecast VT 77111000 
convective heating (°C/day) before 
changes to model constants (see text) 
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Figure lOd. Typhoon KIM 12-h forecast VT 77111000 
convective heating (°C/day) using 
V-bogus 1 and Rh-bogus 1 after changes 
to model constants (see text) 
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Figure 11. Typhoon KIM 24-h forecast VT 77111012 
SLP-1000 mb and winds (m s“^) at 
a = .95 without convective heating effects 
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Figure 12a. 



Typhoon KIM initial SLE-1000 mb and 
winds (m s“^) at a = .95 using 
V-bogus 1 
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Figure 12b. Typhoon KIM initial SLP-1000 mb and 
winds (m s"^) at a = .95 using 
V-bogus 3 
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Figure 12c. Typhoon KIM initial winds (m s~^) at 
400 mb using V-bogus 1, and relative 
humidity (Rh-bogus 1) 



62 




Figure 12d. 



Typhoon KIM initial winds (in s ) 
at 400 mb using V-bogus 3 




Figure 13. Typhoon KIM 12-h forecast VT 77111000 
SLP-1000 mb and winds at a = .95 using 
V-bogus 3 and Rh-bogus 2 
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Figure 14. Typhoon KIM 12-h forecast VT 77111000 
accumulated convective precipitation 
(cm) using V-bogus 3 and Rh-bogus 2 
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Figure 15. Typhoon KIM 12-h forecast VT 77111000 
convective heating C’C/day) using 
V-bogus 3 and Rh-bogus 2 
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Fig, 16. Forecast tracks of Typlioon KIM( 7?1 1 091 2) 

ftpr V-bogus 1 vs V“bogus 3 using Rh-bogus 2 
in the Penn State model. 
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Longitude 

Fig. 17 . Forecast tracks of Typhoon 0LIVE(?80421 1 2) 
using the Penn State(PSU), split semi~impl 
(SSI), and Channel(Ch) models. 
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